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ABSTRACT 

We consider the effects of an outflow on radiation escaping from the infalling envelope around a 
massive protostar. Using numerical radiative transfer calculations, we show that outflows with prop- 
erties comparable to those observed around massive stars lead to significant anisotropy in the stellar 
radiation field, which greatly reduces the radiation pressure experienced by gas in the infalling enve- 
lope. This means that radiation pressure is a much less significant barrier to massive star formation 
than has previously been thought. 

Subject headings: accretion, accretion disks — radiative transfer — stars: formation — stars: winds, 
outflows 



1. INTRODUCTION 

Stars with masses ^ 2{)Mq have short Kelvin times 
that enable them t o reach the main sequence whil e they 

' are still accreting (|Shu. Adams, and Lizanolll987f) . The 
resulting nuclear burning leads to a huge luminosity, 
which produces a correspondingly large radiation pres- 
sure force on dust grains suspended in the incoming gas. 
This force can exceed the star's gravitational pull, pos- 
sibly halting accretion and setting an upper limit on 
the star's final mass. Early spherically symmetric cal- 
culations suggested that this phenomenon s ets an upper 
limit on stellar masses of ~ 20 - 40Mq, llKah3ll97l 

, iWolfire and CassinelHl987l) for typical galactic metallic- 
ities. More recent non-spherical calculations have loos- 
ened t hat constraint by considering the role of accretion 

, disks fNakano 1989; Nakano. Has eeawa. and NormanI 
IT995: .liji.na_^and A dams 19.96). Disks reduce the effects 
of radiation pressure by concentrating the incoming mat- 
ter into a smaller solid angle, increasing its ram pres- 
sure. They also absorb stellar radiation in a thin ring 
and re-radiate it isotropically, casting a shadow of re- 
duced radiation pressure. Even with a disk, however, 
radiatio n pressure can still b e a s ignificant barrier to ac- 
cretion. I.Tiiina and Adarrisl l) 199(11 were able to form a 
50 Mq star only if it had a thin accretion disk with a 
radius ~ 4000 AU. Simulations bv lYorke and Sonnhalteil 
(p002) in 2D found hmiting masses of ~ 40Mq before 
radiation pressure reversed the inflow. However, obser- 
vations show that considerably more massive stars exist, 
and their formation mechanism remains uncertain. 

Recent observations of massive protostars have 
added a n element to this picture. Bcu ther et all 
l|2002albD . iBeuther. Schilke. and Sta^ |2603(l and 
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IBeuther. Schilke. and GuethI 1)20041) report interferomet- 
ric measurements showing outflows from massive proto- 
stars with collimation factors of ^ 2 up to ~ 10. They 
conclude that high and low mass outflo ws have similar 
collimation factors, typically 2 — 5 (Bachil i€ill99(it) . (The 
collima tion factor is the rat io of the outflow's length to its 
width.) iRicher et all l)2000D show that the momentum of 
CO outflows driven by massive stars scales with the bolo- 
metric luminosity of the source in the same manner as 
for low mass stars. From these observations, the natural 
conclusion is that low and high mass stellar outflows have 
a common driving mechanism and similar morphologies. 

Previous theoretical work on massive star formation 
did not include the effects of outflows, and therefore as- 
sumed that the stellar radiation was either isotropic (e.g. 
Jijina and Adams 1996) or had only those anisotropics in- 
duced by the presence of a disk (e.g. Yorke and Sonnhal- 
tcr 2002). In this Letter, we calculate radiative transfer 
through dense envelopes accreting onto massive proto- 
stars, and we study the effects of outflow cavities in the 
envelopes on the radiation fleld. In § |21 we present our 
models for radiative transfer and for the protostellar en- 
vironment. We then present the results of our calcula- 
tions in § 121 and discuss conclusions in § 0] 

2. MODELS 

2.1. Radiative Transfer Methodology 

We used the Monte Carlo /diffusion r adiative transfer 
code written bv lWhitnev et all l)2003albfl to find the tem- 
perature distribution of gas in a circumstellar envelope. 
The opacity comes from dust grains, which are ther- 
mally well-coupled to the gas in the high density envelope 
(ISpitzer 1978). The grain size distribution depends on 
the location in the flow. The stellar spectrum is a Ku- 
rucz model atmosphere. Rather than discuss the code, 
grain, and stellar spectral models here, we refer readers 
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to lWhitnev et all l|2003aj) . and to 5 12 .31 for a discussion of 
opacity in the outflow.We used 10^ photons, sufficient to 
produce high signal-to-noise spectra. 

To calculate the radiation force on the gas, we be- 
gin with the temperature distribution determined by 
the Monte Carlo code. From a given position, we 
pick ray directions based on the HEALPix scheme 
()G6rski. Hivon. fc Wandel^ Il999|) . Along each ray, we 
solve the transfer equation for emission and absorption 
by dust grains to obtain 1^ (n) , the specific intensity com- 
ing from direction n. We also perform this calculation 
along a ray to the star itself. Once we have found the 
intensity along each ray, we compute the flux 
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where Arrays is the number of rays. We increase Arrays 
until the change in flux is less than 2.5% between iter- 
ations. We then integrate over frequency to obtain the 
radiation force per unit mass, frad = / k,^F^ dv. 

We ignore scattering when calculating the force (but 
not the temperature). Scattering of IR photons is neg- 
ligible. UV and visible photons can scatter significantly, 
but for the fiducial envelope and star we use, photons at 
the peak of the stellar spectrum are all reprocessed into 
IR in the inner few AU of the envelope. (In the disk the 
distance is vastly smaller.) Outside this layer, neglecting 
scattering does not change the radiation pressure force. 
For the same reason, the choice of stellar spectrum model 
is relevant only in the inner few AU of the envelope. 

2.2. Star, Disk, and Envelope Properties 

To study the effects of outflow cavities, we choose a 
single fiducial model for the star, disk, and envelope, and 
vary the properties of the outflow. We place an M, = 50 
Mq zero-age main s equence star in a 50 Mq envelope. 
The ZAMS models of lTout et aill|1996|) . predict a radius, 
surface temperature, and luminosity o f fij. = 10.8 
T; = 4 .3xlO-^K,andL, = 3.5 xlO^ L^. iMcKee and TanI 
()2003j) predict a formation time of ~ 10^ yr for such a 
star, so we adopt an accretion rate of M^, = 5 x 10"'* AIq 
yr~^. The accretion luminosity is Lace ~ GK'UM^:/ R^, — 
7.1 X lO^i©, negligible in comparison to the central star. 

Although the cores that form mass ive stars are prob- 
ably turbulent ijMcKee and Tan 2003), for simplicity we 
adopt a si mple rotationally flattened density distributio n 
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Here P2 is the Legendre polynomial, r and 9 give the po- 
sition in the envelope in spherical coordinates, and i?cen 
is the centrifugal radius of the flow. The gas at position 
(r, 6) was at angle 6*0 when it began falling toward the 
star, where 9o is given implicitly by 10} . The centrifu- 
gal radius i?ccn is determined by the angular momentum 



of material arriving at the star. When the last material 
in the core accretes, it is roughly i?ccn — PRcore, where 
P is the ratio of the core's rotational kinetic energy to 

C ravitational binding e nergy and for most cores /3 ~ 0.02 
Goodman et a]llT993) . Since our core is half accreted, 
we take i?cen ~ PRcoie/2. We therefore adopt values of 
^coro = 0.22 pc and i?cen = 400 AU, which give the cor- 
rect mass in the envelope and satisfy i?ccn ~ /9^coro/2. 

In addition to the envelope, our star has a disk of ra- 
dius i?cGn- Observations of disks around massive stars 
are limited, and we therefore adopt a disk mass of 
Md = M*/10 — 5Mq based on the consideration that 
a disk more massive of order 10% of the central ob- 
j ect is likely to be subject to gravitational instabilities 
ijShu et all ll9901 that will cause matter to accrete until 
marginal stability is restored. We model the disk sur- 
face density and scale height as power laws E oc r~^-^^, 
h = 0.035 {r/R^) AU. We base the choice /i ex r on 
images from simulations of massive star formation (e.g. 
Yorke & Sonnhalter 2002) which show disks with little 
or no flaring. 

2.3. Outflow Properties 

We can partially describe the outflow cavity from ob- 
servations. As discussed in § ^ high mass outflows have 
collimation factors from 2 — 5, similar to low mass stars. 
This corresponds to half-angles from — 15°, including 
the uncertain inclination. We therefore test opening an- 
gles of 9o = 5°, 10° and 15°. Unfortunately, observa- 
tions constrain only the asymptotic opening angle when 
the outflow is far from its parent core. Cavity walls 
from nearby low mass sources seen in near-IR are gen- 
erally curved when viewed on sufRciently small scales 
(e.g. Padgett et al 1999 and Reipurth et al 2000). This 
curvature is a natural result of out flow collimation by 
the rotationally-flattened envelope (Wilkin and Stable^ 
^998, 2 003) . For massive protostars that are still embed- 
ded in dense envelopes, near-IR observations are impossi- 
ble due to dust extinction. Millimeter observations, even 
with interferometers, are unable to probe length scales 
comparable to the size of the protostellar disk, at which 
we exp ect the stron gest collimation and curvature. Fol- 
lowing [Whitne^^Fal (2003b), therefore, we parameter- 
ize the uncertain shape of the cavity as z = aw^ , where 
z is the vertical distance from the star, 137 is the distance 
from the outflow axis, and a — i?corc cos 9o/ sin'' 9o is a 
constant chosen to give an opening angle of 9o at the 
edge of the core. We try values of 6 = 1.25, 1.5, and 
2.0, coupled with a fixed wind opening angle 9o = 10°, 
to study the effects of variations in cavity shape. As a 
baseline, we also study a case with no outflow cavity. 

Some gas at the outflow base will be ionized by the 
star's UV flux (Tan & McKee, in preparation). Whether 
the ionized gas remains near the star depends on the out- 
flow's structure at its base, where the density is highest 
and recombinations are fastest. If the outflow is ionized, 
it will be too hot to contain dust grains, and its opacity 
will come primarily from resonant scattering by metal 
ions (jCastor. Abbott, fc Kleinlll975|) . This is a complex 
subject beyond the scope of this Letter. We simply note 
that the opacity produced by this scattering must be far 
smaller than that produced by dust grains in the enve- 
lope. In the case that the outflow cavity is ionized, we 
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Fig. 1. — The color maps show the gas temperatures for each 
of our models. The models are (a) no wind cavity, (b) 8o =5°, 
b = 1.5, (c) Oo = 10°, b = 1.5, (d) Bo = 15°, b = 1.5, (e) Oo = 10°, 
b = 1.25, and (f) = 10°, 6 = 2.0. The red dots inside the cavity 
in panel (c) are the result of a minor code bug. 



may therefore set the opacity in the cavity to zero. 

If the outflow is neutral, grains will re-form, growing 
in radius r^^ at rate (cf. Hoyle 1946) 

Tgr = avtp/{4,pgr), (5) 

where a is the mass fraction of the element of which the 
grain is composed, Vt is the thermal velocity of gaseous 
atoms of that element, p is the gas density, and is 
the grain density. Growth is fastest for carbon grains, 
since carbon is the most abundant refractory metal in 
the galactic ISM an d has a large thermal ve locity due to 
its low atomic mass. iSofia and Meved 1)200 If) estimate its 
mass fraction to be ac « 3 x 10~^. Grains can only con- 
dense below the dust destruction temperature of Tdust ~ 
1600 K, so vt'^l km s^^. The density of carbon grains 
is ~ 1 g cm~^. For our adopted cavity shape, the wind 
density at a distance z from the equatorial plane must 
be roughly fu)M^/[2TTVw (z/a)^^''], where Vw is the wind 
velocity (taken to be constant) and fw is the fraction of 
mass reaching the star that is ejected into the wind. We 
take /„ = 0.2 and = vk = {GM^/R/)^/^ « 1000 
km s-i (|Richer et all 120001) . where vk is the Keplerian 
velocity at the stellar surface. 

The smallest distance at which grains can form is 
the dust destruction radius, which we calculate to be 
-Rdust ~ 48 AU for our fiducial model. We may therefore 
integrate lO from i?dust to i?corc to obtain the maximum 
size that grains can attain before escaping the core. The 
largest grain sizes occur for 9 — b° and b — 1.5, which 
gives Tgr < 1.1 X 10~^ /im - in effect, grains cannot grow 
past the size of molecules. In contrast, the typical grain 
sizes in the envelope are 0.1 — 0.2 /im. Stellar and enve- 
lope radiation has wavelengths A ^ 27rrgr, so grain opac- 
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Fig. 2. — The plots show the force per unit mass due to radiation 
and gravity {black line), at angles of 6 = 60°, 75°, and 90° from 
the polar axis. The plots begin at the minimum value of r that is 
inside the infalling envelope. The three left panels show radiation 
with no cavity (red line), do = 5° (green line), do = 10° (blue 
line), and do = 15° (purple line), all with b = 1.5. The three right 
panels show radiation with no cavity, b = 1.25 (green line), b = 1.5 
(blue line), and b = 2.0 (purple line), all with do = 10°. 



ity scales as k oc r^^. Thus, even if all the C goes into 
grains, the opacity in the outflow is smaller than that 
in the envelope by at least a factor of ^ 800. If we in- 
cluded cooling of the outflow and evaluated the grain size 
at radii smaller than i?corc, the reduction would be even 
larger. We conclude, therefore, that dust opacity in a 
neutral outflow is negligible in comparison to the opac- 
ity in the infalling envelope. Since for either a neutral 
or an ionized outflow the opacity is negligible, we set the 
opacity in the outflow cavity to zero in our calculations. 

3. RESULTS AND DISCUSSION 

We show the gas temperature distribution for each of 
our models in Figure ^ Although our core extends to 
0.22 pc, we concentrate on the inner few thousand AU, 
where the radiation pressure force is strongest. The plot 
shows that, in models with a cavity, the gas is hotter 
on the inside edge of the outflow cavity but is cooler 
elsewhere. Since, due to rotational flattening, most of 
the gas is attempting to accrete close to the equatorial 
plane, this means that the bulk of the accreting gas is 
cooler in runs with an outflow. The effect is stronger 
for outflow cavities with larger opening angles, and for 
cavities that are wider at their base. 

We show the radiation pressure force versus radius at 
angles of 6* = 60°, 75°, and 90° from the pole in Figure 
121 Consistent with the reduced temperature, the radia- 
tion pressure force is also substantially smaller in runs 
with an outflow cavity. For our "intermediate" model 
with 6o = 10°, b = 1.5, the radiation pressure force is 
smaller than in the run with no cavity by a factor of up 
to 4.6. For the other runs, the peak reductions in radi- 
ation force were factor of 1.7 {do = 5° , b — 1.5), 14.4 
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[Oo = 15°, b = 1.5), 4.2 {Oo = 10°, b = 1.25), and 7.2 
(0o = 10°, 6 = 2.0). Thus, the amount by which the wind 
cavity reduces radiation pressure shows moderate depen- 
dence on the wind curvature and opening angle, but for 
only one combination of parameters was the reduction 
less than a factor of 4. Comparing radiation and gravi- 
tational forces, it is clear that this reduction can mean 
the difference between accretion halting or continuing. 

To check how the results depend on properties of the 
envelope, which has a comparatively low surface den- 
sity compared to most high-mass cores, we re-ran the 
intermediate and no cavity cases with identical i?cen and 
i?coro, but a 100 Mq envelope. The case with no wind 
cavity showed little change in force with envelope mass. 
With a cavity, the increased envelope mass decreased the 
radiation pressure force at equatorial angles so that radi- 
ation was stronger than gravity only aX 9 1^ 60°, increas- 
ing the fraction of solid angle through which accretion 
could occur. Thus, our results likely represent a lower 
limit on the reduction in radiation pressure force that 
outflow cavities actually produce. 

To determine how the effect of the cavity compares to 
that produced by the disk and the rotationally-flattened 
envelope alone, we also considered a spherical envelope, 
with no cavity, no disk, and -Keen = in l(2l-l@I- We 
compared this to the run with our fiducial disk and en- 
velope parameters and no cavity. At angles 1^ 75°, the 
radiation pressure force in the spherical case was a fac- 
tor of a few higher, while at angles ^ 60° it was com- 
parable or smaller. Thus, our intermediate wind cavity 
reduces the radiation pressure force relative to the disk- 
and-envelope only case by about the same amount that 
the disk-and-envclope case reduces it relative to purely 
spherical. CoUimation of the radiation field by the disk 
and envelope and coUimation by the outflow cavity are 
about equally important, and reinforce each other. To- 
gether, they reduce the radiation pressure force in the 
intermediate cavity case by a factor of 10 relative what 
one would flnd for an isotropic radiation field. 

In all our tests, the degree of radiation coUimation is 
roughly consistent with expectations. Tan & McKee (in 
preparation) show that the fraction of radiation escap- 
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ing an envelope through a path of optical density r is 
proportional to (1 -f t)~^. In our case, the flux fraction 
escaping from the core at angles between 6 and 6 + d9 
should roughly satisfy 

dF{9) o^[l + T{9)r^d[cos{9)], (6) 

where t{9) is the optical depth at angle 9. There is con- 
siderably uncertainty in applying this to our problem be- 
cause t{9) is frequency-dependent. However, if we use 
the Rosseland mean opacity at Tdust = 1600 K, then at 
most angles (O predicts roughly the correct flux fraction. 

4. CONCLUSION 

We have shown that an outflow can substantially 
change the radiation fleld, and radiation pressure, around 
a massive protostar. The outflow cavity provides an op- 
tically thin channel through which radiation can escape, 
significantly reducing the radiation pressure. With no 
wind cavity, in our fiducial model the radiation pressure 
force is stronger than gravity essentially everywhere ex- 
cept inside the accretion disk, and it is therefore likely 
that accretion would be halted. In our intermediate 
model, outside the centrifugal radius radiation is weaker 
than gravity over about tt sr, providing a large funnel 
through which accretion can continue. The calculation 
we have performed here is only a proof of principle, and 
we are currently performing 3-D radiation hydrodynamic 
AMR simultions. Our results strongly argue, however, 
that the presence of outfiows provides a mechanism for 
circumventing the radiation pressure limit to protostel- 
lar accretion. Surprisingly, outfiows that drive gas out of 
a collapsing envelope may increase rather than decrease 
the size of the final, massive star. 
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